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ABSTRACT 

We present a method to implement the idea of iJain fc Tavloil l)2003ft to constrain cosmological pa- 
rameters with weak gravitational lensing. Photometric redshift information on foreground galaxies is 
used to produce templates of the mass structure at foreground slices zg, and the predicted distortion 
field is cross-correlated with the measured shapes of sources at redshift Zg . The variation of the cross- 
correlation with Zs depends purely on ratios of angular diameter distances. We propose a formalism 
for such an analysis that makes use of all foreground-background redshift pairs, and derive the Fisher 
uncertainties on the dark energy parameters that would result from such a survey. Surveys from the 
proposed SNAP satellite or the LSST observatory could constrain the dark energy equation of state 
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0.035/g]^y after application of a practical prior on f2„j. Advantages of 
this method over power-spectrum measurements are that it is unaffected by residual PSF distortions, 
is not limited by sample- variance, and can use non-linear mass structures to constrain cosmology. The 
signal is, however, very small, amounting to a change of a few parts in 10'^ of the lensing distortion. 
In order to realize the full sensitivity to cosmological parameters, the calibration of lensing distortion 
must be independent of redshift to comparable levels, and photometric redshifts must be similarly 
free of bias. Both of these tasks require substantial advance over the present state of the art, but we 
discuss how such accurate calibrations might be achieved using internal consistency tests. Elimination 
of redshift bias would require spectroscopic redshifts of ~ 10'' — 10^ high redshift galaxies — fewer for 
lensing surveys less ambitious than SNAP or LSST. 
Subject headings: gravitational lensing; cosmological parameters 



1. INTRODUCTION 

Weak gravitational lensing is already one of the 
more accurate constraints on cosmological parameters: 
several groups have measured the power spectrum of 
the shear induced on background galaxies by fore- 
ground mass fluctuations, leading to « 10% con- 
straints on the matter power spectrum normalization 
(T8 (Bacon et al. 2002; Brown et al. 2003; Haman a et al. 
I2p02; Hoekstra, Yec. & Gladders 200 2t I.Tarvis et al.\ 
12003.: Befregier et al. .2002: .van Waerbeke et a.l\ I2nn2|) . 
The weak lensing method has potential for much more 
precise constraints on trg and other parameter combi- 
nations, once larger sky areas are surveyed. The addi- 
tion of redshift information on the source galaxies allows 
measurement of the mass power spectrum as a func- 
tion of redshift, which is expected to greatly increase 
the accuracy of weak lensing cosmological constraints 
(lHulir99 9. 2002a b; Huterer 2002; Abazaiian & Dodelsa 
120021 lHeaven£ll2003t IRefreeier et ali I2003t IKnoxll20d" 
and permit r econstruction of the three-dimensio nal mass 
distribution ijTavlod l200lt iHu fc KeetonI 1200^ . Weak 
lensing surveys currently underway such as the CFHT 
Legacy Survey^ (CFHLS) and the Deep Lens Survey^ 
will gather photometric redshift information to facilitate 
this "tomographic" analysis. 

The successful constraint of the mass power spectrum 
and cosmological parameters at « 1% levels will require 
reduction of several sources of systematic error. On the 
theoretical side, there are not yet predictions of the mass 
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power spectrum P{k,z) that are acc urate to the per- 
cent level in the non- linear regime fsee lLinder fc Jenkinj 
(1^003) for further discussion). A concerted application 
of iV-body computing would likely yield the dark-matter 
spectrum to desired accuracy, though on very small scales 
the contribution of baryons to the power spectrum must 
be included, and will be difficult to calculate. The inac- 
curacies of power-spectrum estimation may be bypassed 
by using only large-scale information, where linear or per- 
turbative calculations suffice, but this means discarding 
most of the lensing information, which lies at non-linear 
scales. 

On the measurement side, currently published power- 
spectrum measurements all show contamination by sys- 
tematic errors at the « 10% level or higher. This 
systematic power is likely residual from the process 
of correcting galaxy shapes for point-spread function 
PSF) effects. New metho dologies ha ve been introduced 
Bernstein fc Jarvisll200l iRefregied l2003t iKaised l2000j) 
which should greatly reduce the systematic contamina- 
tion, but this remains to be demonstrated. There are also 
subtle difficult ies in calibrating the weak lensing signal to 
1% a ccuracy l|Bernstein fc Jarvisll200l iHirata fc SeliakI 



iJain fc TavTo^ l|2003D introduce a new method for 
analysing weak-lensing data with depth information 
which promises to largely bypass these systematic dif- 
ficulties. The basic concept is to use the survey's pho- 
tometric redshift data to create a map of the foreground 
galaxies, from which an estimated map of the foreground 
mass can be made. This foreground mass slice induces 
shear on all the galaxies in the background. The mea- 
surement consists of tracking the amplitude of the in- 
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duced shear as a function of the background redshift. 
Once a foreground shear "template" at zi is created, the 
dependence of shear on the source redshift Zg is given 
solely by geometric factors. By taking ratios of induced 
shears at different z^ , any sensitivity to errors in the fore- 
ground mass template is cancelled out, leaving us with 
a purely geometric observable. One is not attempting to 
discern the line-of-sight structure, so this is not tomogra- 
phy. We could call this "cross-correlation cosmography" 
since the measured quantities are metric elements of the 
homogeneous cosmology. 

The idea of using lensing effects on sources at 
differing redshift to constrain cosmological parame- 
ters has been examined befo re jLi nk & Pierce 1998; 
iGautret. Fort, fc Mellieil l2000t [Colse. Kneib. fc Soucaill 
12002: ,SerenQ 2QQ3). with Jain fc Taylor offering a 
methodology for massive datasets that is insensitive to 
the nature of cluster or halo profiles and examines a time- 
dependent dark energy component. 

The practical advantages of cross-correlation cosmog- 
raphy include: 

• No theoretical estimate of the mass power spec- 
trum is required. The mass fluctuations are esti- 
mated directly from the foreground galaxy distri- 
bution, and any inaccuracies {e.g. bias) are can- 
celled out in the analysis. The observables are 
hence calculable to arbitrary accuracy. 

• With no need for a modelled power spectrum, we 
can use shear information on all scales sampled by 
the data. The lensing distortion variance {(P) is 
~ (5%)^, which is at least two orders of magnitude 
larger than the variance due solely to linear-regime 
fluctuations. Hence the signal-to-noise ratio [S/N] 
of the weak lensing data will be much higher. 

• Because the background shear is being cross- 
correlated with a template, the shear signal en- 
ters only linearly into the statistics. Systematic 
power {e.g. PSF residuals) average to zero since 
they will not correlate with the foreground shear 
template. This makes the PSF-correction task im- 
mensely easier than for power-spectrum tomogra- 
phy, in which systematic power adds to the shear 
power. For cross-correlation cosmography, system- 
atic power may increase uncertainties but does not 
bias the results. In the language of experimenters, 
we have changed from a total-power measurement 
to a phase-sensitive method. 

• Uncertainties in calibration of the shear signal also 
cancel if they are independent of Zg, since we will 
be interested only in ratios of shear at different Zg. 
We shall see below, however, that the method is 
highly sensitive to differential calibration errors. 

Taking the ratio of shear signals at different Zg provides a 
pure geometric measurement, but taking this ratio does 
decrease the dependence of the signal on cosmological 
parameters. Large numbers of galaxies must be surveyed 
in order to reduce random errors to make up for the 
smaller signals, but surveys of up to 10^ galaxies are 
currently being planned. 

In this paper we investigate the potential of cross- 
correlation cosmography by proposing an analysis 



methodolgy that appears close to optimal, and calcu- 
lating the expected cosmological-parameter uncertainties 
with the Fisher-matrix method. We then apply this anal- 
ysis to some planned weak lensing surveys. Finally we 
return to the issue of systematic errors, exploring the 
accuracies that will be required for photometric redshift 
estimates and shear calibration. 

2. CROSS-CORRELATION FORMALISM 

2.1. The Observable Quantity 

The underlying assumption of the method is that 
the observed ellipticity of background galaxy k £ 
{1,2,. ..,7V} is 

ef"" = ^ de{xk)gek + dk,sys + efc, (1) 
The quantities of interest will be the geometric factors 



where Xi and Xs are the comoving distances to a lens and 
source planes, and r{x) is comoving angular diameter dis- 
tance. For a flat Universe, r{x) = X- Each of the is 
the distortion field'^ imparted on a source plane at x = c» 
by the mass in redshift shell £. It can be expressed as an 
integral over the mass distrib ution in the red shift shell 
by we ll-known formulae, e.g. IBartelmann fc Schneide'3 
l)2001t) . Note that the distortion d is twice the shear 7 
in the weak- lensing limit. Its variance (cfj) can be calcu- 
lated as an integral of the mass power spectrum P{k, zi), 
with the formula in Appendix^ Our method, however, 
will be to make as little use as possible of the power 
spectrum. 

Equation is correct only in the limit that the in- 
duced distortion df, the instrumental systematic distor- 
tion dfe sys, and the intrinsic galaxy shape are all <C 1. 
The ad dition operator in distort ion space is in fact non- 
linear ijBernstein fc Jary is" 2002), a detail that will be 
important, but manageable if we avoid the « 1% of the 
sky where the lensing distortion is strong. 

Any weak lensing measurement maximizes the likeli- 
hood of the measured e^'''' under some model for the 
shear planes d^ (and perhaps the systematic error con- 
tribution simultaneously). In power-spectrum tomogra- 
phy, one varies the assumed power spectra of the d^ to 
best match the observed covariances of the e^.'^''. In our 
case, we will assume that is possible to produce some 
estimate d^ of the distortion pattern without reference 
to the background galaxy-shape information. In prac- 
tice this template distortion would be developed from the 
galaxies identified near zt from the photometric redshift 
data. This could be as simple as convolving the galaxy 
distribution with an isothermal halo model, or perhaps 
involve a more sophisticated identification of groups and 
clusters. We allow that the template distortion is inac- 
curate to some level, due to galaxy bias or other errors: 

(d, -d,) =/3,(d^) (3) 

is a measure of the fidelity of our template, and is 
related to the bias b and correlation coefRcient r between 

^ The distortion denoted as 6 in IBernstein &: Jarvi IS 
labelled d here in order to avoid confusion with the Kronecker 
delta. 
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galaxie s and mass, e.g. as measured bv iHoekstra et aD 

To isolate the geometric term, we correlate the tem- 
plate distortion with the measured distortion. We wish 
to sum the measurements from all source galaxies in some 
bin s of redshift centered on Zs . When the shape noise is 
dominant, the sum which is optimal in the sense of best 
S/N on gis is 



.obs 



(4) 



fees 



Here Ns is the number of source galaxies in bin s. Under 
the model in Equation the X^s will be 



Rff, EE (df • Af,) 



(d, • d 



sys/s 



(5) 
(6) 



Here the subscript s on the average means that we are 
averaging over the galaxies in source bin s, as in equa- 
tion^ not over realizations of the distortion fields. For a 
given survey, the cosmology predicts the gi^ values, the 
templates are known, and the R^i are free parameters. 
The cosmological parameters are determined by fit- 
ting the X — {Xiis\ data vector to the model © with 
free parameters {pi,Rui}, then marginalizing over the 
cross-correlations {Ru'} to obtain confidence bounds for 
the Pi . The intrinsic ellipticities and perhaps systematic 
errors act as measurement noise on each Xgg . 

We now assume that the intrinsic ellipticities and 
the systematic-error distortion d^^sys are independent of 
the template. We also assume that the distortions and 
templates are uncorrelated between lens shells, so that 
{Ru>) = for £ ^ This last assumption is true if the 
shells are much thicker than the correlation length of the 
mass distribution; any correlations could be accommo- 
dated in a more sophisticated analysis, or eliminated by 
applying high-pass filters to the templates. With these 
assumptions, the expectation value for the Xis when we 
average over realizations of the distortion fields is 



(Xis) = {Ru)gts = Pi{dl)gis 



2.2. Covariance Matrix and Likelihood 



(7) 



In order to calculate probabilities for fits to the X vec- 
tor we need its covariance matrix C^^. From the definition 
Q and the model we can obtain the covariance in a 
straightforward manner. Note that when fitting to the 
model Q we will be considering only the shape noise 
and systematic errors to be random variables. The dis- 
tortion planes, or more precisely the cross-correlations 
are considered as free parameters in the fit, not ran- 
dom variables. We will henceforth ignore the systematic 
distortion dgys, assuming that it is a minor contributor 
to the noise compared to the shape noise e. We have 
already seen that the systematic errors do not bias the 
observables if they are uncorrelated with the templates. 

We assume the to be independent, with Var(e+) = 
Var(ex) = the variance in each component. Intrin- 
sic correlations of galaxy shapes may slightly inflate our 
uncertainties, but should again be uncorrelated with the 



dp and hence will not bias the fit. In this case we have 



£ = f' 

' i^e 



(8) 



The I' cross-correlation is negligibly small if the sur- 
vey is much larger than the correlation length of the dis- 
tortion templates. 

The distributions of the d^ and are very non- 
Gaussian. For a survey with many source galaxies and an 
area much larger than the distortion correlation length, 
however, the distribution for Xg^ should approach Gaus- 
sian by the central limit theorem. The likelihood of an 
observed X vector given a choice of model parameters is 
then given by the usual Gaussian formula with covari- 
ance matrix C^;. With the likelihood, we may produce 
confidence bounds on the parameters of the underlying 
geometry, once we marginalize over the unknown . 

2.3. Fisher Uncertainties 

Given the likelihood function, the minimal uncertain- 
ties on the model parameters may be derived from the 
Fisher matrix: 



F — — 



dpidpj 



(9) 



where C = logL. The Fisher matrix for the Gaus- 
sian ca se is well known, e.o. lTegmark. Tavlor. & HeavensI 

F,,. iTr(C-iC,,C-iC,,-) + (X;^>C-i(X,,.). (10) 

As usual the subscripts with commas denote derivatives. 
This simplifies considerably in our case, as C^: is diagonal 
and independent of the parameters. 

The model parameters are of two kinds: the nuisance 
parameters {Rw} and the cosmological parameters {pi\. 
We wish to split the Fisher matrix into submatrices V rr, 
FpR, and Fpp. Henceforth we will adopt the convention 
that Greek indices (or i) range over the lens planes, the 
s index runs over source planes, and other Latin indices 
(except i) range over the cosmological parameters. Keep- 
ing in mind that it takes two Greek indices to specify one 
R component, the three submatrices Frr, FpR, and Fpp 
are respectively 



F aPfi'y — 
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-gpsg-is 



N 



olkdl 
N 



-G 



Pi 



' iaf3 ' 



<^i{dl 



y^i?a-y(Gi)/3^ 



N 



72 \ ^ 



7 ^ RapRaiiGij)^^. 



(11) 

(12) 
(13) 
(14) 



The geometric factors are encoded in the matrices 

dzs-^gasgps (15) 
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{Gi)al3 = 



dzs dpi 



dn dgas dgps 
' dzs dpi dpj 



(16) 
(17) 



The redshift distribution of sources dn/ dz is normalized 
to unit integral, and TV is the total number of sources in 
the survey. There is no loss of information in moving to 
infinitesimal source-redshift bins. 

2.4. Prior Information 

We are interested in the covariance matrix Cp of the 
geometric parameters after marginalization over the R 
parameters. If L is the number of lens planes in the 
model, there are free Rw being marginalized, and 
after this marginalization the constraints on the cosmo- 
logical parameters are relatively weak for envisioned sur- 
veys. 

There is, however, additional information about the 
R values that has not been incorporated into the model 
of Equation JSJ, and hence not in the Fisher matrix of 
Equation (|lll) . We expect that the Rw are near zero 
for I ^ i' because the lens planes are uncorrelated. We 
incorporate this knowledge with a prior probability on 
the R values. If the prior is Gaussian, then we may 
add the inverse of the R covariance matrix to the Yrr 
component of the Fisher matrix. 

Note that until this point we have had no use for an en- 
semble of mass distributions; all of the quantities in the 
solution and Fisher matrix for the cosmological param- 
eters make use of just the Rw for the single realization 
of the mass distribution that exists in our survey field. 
It is only in assembling a prior on these values (and in 
estimating typical values in t|2.5|l that we make use of 
the ensemble properties of the mass distribution. This is 



why the sample variance contribution is unimportant in 
the uncertainties of the cross-correlation method. 

In the Appendix we show that the covariance matrix 
for Rill is nearly diagonal. The diagonal terms are 



Q/3 



(18) 



Here is the total solid area of the survey and Q.af} is 
a measure of the area of coherence of the two distortion 
fields da and dp. For redshift ranges of interest, fia/J ^ 
1 arcmin^ . The only off-diagonal terms in the covariance 
matrix for the R values are Cov(i?jj^i?^jj). To simplify 
the following algebra, we take the covariance matrix for 
our prior on R to have diagonal elements that are twice as 
large as Equation (|18|l and drop the covariance between 
Raf} and Rpa- This prior distribution is less restrictive 
than is the full covariance matrix, so we are at liberty to 
make this choice. 

As a cautionary step we do not make use of the prior 
information on Raa- This is because Raa appears in 
the expectation value {X^s) in Equation ([7|), and hence 
prior assumptions may bias the fit for the cosmological 
parameters that drive the gas in this fit. 

With this simplified, diagonal prior for the R values, 
the Fisher matrix for the system is altered from Equa- 
tion l|ll|l as 



= 5aL 



N 



-I- S/3j{l — Sa/s) 



2nQ.al3{d?ij) 



{G + Pa. 



(3-1 



(19) 



Pq is simple, having only diagonal elements and a zero 
at the a element of the diagonal. Here n = N/Q, is the 
sky density of source galaxies. 



J 



2.5. Marginalization 

We now return to the marginalization over the Rap to obtain a covariance matrix Cp for the cosmological parameters. 
Using a common matrix identity, we obtain 



{Gp ^)y - (Fpp - FpijF^^jFj^) - -J X! 71^7 X! ^<^P^<^i i^v ~ (G -I- P^) ^G 



(20) 



This gives the Fisher uncertainties on cosmological parameters for a survey over a sky with given distortion field and 
hence given R values. We next average over realizations of the mass distribution in the Universe, which requires that 
we calculate 



(RafsRaj) = Gov{Ra0Raj) + {Rap) {Raj) = ^/S-) 



^2\/j2\ , r /n \2 
r Oap{-n.aa) 



{di){di)- 



(21) 



Combining Equations (0), H2Q(I . and (|2HI . we obtain the expected (inverse) covariance matrix for the cosmological 
parameters: 



(c;^) 



N 



f3l{dl)(G.,-Gj{G + Pa)''G, 



E(4> 

/3 



Q/3 



n 



G, 



Gf(G 



'/3/3 



(22) 



We note first that there is no "sample variance limit" to 
the cross-correlation cosmography measure of the cosmo- 
logical parameters — uncertainties scale as 1/\/N as the 
density of sources becomes large. Sample variance typi- 



cally arises when one does not measure enough indepen- 
dent patches of the density field to adequately constrain 
its power spectrum. In our application, however, we are 
not using the power spectrum to measure cosmology, we 
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instead are just cross-correlating the observations with 
whatever mass distribution the Universe gives us. Hence 
there is no sample variance contribution. 

The left-hand term in Equation (|22|l arises purely from 
the shape noise in the measured field that is being cross- 
correlated with the templates, while the right-hand term 
arises because the distortion from one lens slice acts as 
a correlated noise source for all the Xig. The right- 
hand term has an amplitude « L/N^ relative to the 
left-hand term, where L is the number of lens planes 
and Nc — U,/VLa_i3 is the number of independent cells or 
patches of the distortion fields. In almost any large sur- 
vey we should be able to choose the lens-plane width Az^ 
to be sufficiently narrow that the gi^^ do not vary signifi- 
cantly within the bin, and still have L/Nc <C 1. In other 
words we will usually have more than enough indepen- 
dent patches to decorrelate the different lens planes. 

In the regime where the left-hand term dominates, the 
parameter accuracy is essentially independent of the fore- 
ground binning, as long as the fidelity Pi remains high 
and the shells remain uncorrelated. 

3. APPLICATION TO CANDIDATE SURVEYS 

3.1. Cosmological Parameters 

Any parameter that influences growth factor a(t) or 
curvature of the Universe will have some effect on the 
X data via gis] in particular we will be interested in the 
present-day matter and energy densities ilm and fix , and 
the parameters wq and Wa that a.pproximate t he equation 
of state of the dark energy via llLindeill200r}) 

Px ^[wo+Wa{l~a)]px- (23) 

Wc will henceforth assume a fiat Universe fix = 1 — ^m, 
in which case wc have (Lindcr 2003) 

r[x{z)]=x{z) I dz' [r!„(l + z'f + 

J^o Jo 

(1 - an)(l + z')'(^+""'+"'")e-3"'»^'/(i+^')] 

Unless otherwise noted, we will calculate the sensitivity 
to departures from the canonical ACDM model {flm = 

0. 3, Wo = -1, Wa = 0). 

3.2. Mass Variance 

The distortion power generated by a mass slice at z; 
and the coherence areas flap can be estimated by ap- 
plying the formulae in the Appendix to a model for the 
non-linear evolution of the power spectrum. We reiterate 
that errors in this model will affect only our forecasted 
uncertainty, not the derived cosmological parameters. 

We use the nonlinear m a.ss power spectru m give n 
by the fitting formula of iPeacock fc Dodd^ l)1996D . 
with implementation for weak lensing as described in 

1. Tain fc SeliakI l|1997D . We need this only for the fidu- 
cial ACDM model, for which the differences between dif- 
feren t fits to N-body simulations are small (fSmith et all 
l2002t) . The integrals over the power spectrum are cut off 
for k > 27r/(50kpc); this is a crude way of accounting 
for the fact that we will have to exclude the central re- 
gions of galaxy clusters, where the lensing is strong and 
where the light of the foreground galaxies may preclude 
observation of the background galaxies. The results are 
not sensitive to the choice of the high-A: cutoff. We will 
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Fig. 1. — The top panel shows the fractional change in the 
geometric factor gi^ — {xs — Xe)/Xa when we shift from a 
pure ACDM Universe to one with Wa = 0.2; this should be 
discernible at la significance in the SNAP survey. The hor- 
izontal axis is Zs and each line corresponds to a different Zi. 
The triangle at the end of each line marks ze. Because the 
mass normalization in each lens slice is free to vary, vertical 
shifts of each line carry no cosmological information, so we 
align them all to be unity at zg. The cosmological information 
is carried in the departures of each line from horizontal; these 
departures are small, amounting to only a few parts in 10~^ 
change in the induced background distortion. The smallness 
of this signal implies that the calibration of the distortion 
measurement and the source redshifts must be accurate to 
roughly a part in 10^. The lower panel plots the assumed 
source redshift distribution (dotted line) and the expected 
distortion variance per unit redshift (solid line) using esti- 
mated non-linear power spectra. The dashed line shows the 
relative contribution of different lens planes to the constraint 

on Wa- 



assume that the estimated templates have the same 
variance as the true distortion fields d, but with fidelity 
Pe = 0.8. The uncertainties on cosmological parameters 
will scale as (3^^. 

The bottom panel of Figure ^ plots the strength of 
lensing distortion vs redshift, d{d?)/dz that arises from 
this model. For the coherence angle we use 

n^p = (0.24arcmin)2(x„x/3)"°-', (24) 

which is within ~ 25% of the more carefully calculated 
values. The coherence area has only a minor impact on 
the cosmological uncertainties. 

3.3. Redshift Distribution 

We adopt the common guess for the redshift distribu- 
tion of faint galaxies, 

gocz2exp[-(z/zo)^-^], (25) 

and will select zq and the overall density n to crudely 
mimic the expectations of several future surveys. 

3.4. Parameter Forecasts 

Table ^ lists crudely approximated parameters for 
three possible weak-lensing surveys: the CFHT Legacy 
Survey, just beginning and expected to take 5 years to 
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Fig. 2. — Fisher uncertainty ellipses for dark energy pa- 
rameters derived from three candidate weak lensing surveys 
are plotted. All ellipses are 68% confidence two-dimensional 
regions (A^^ ~ 2.3) after application of a Gaussian prior 
on Qm with a — 0.03 and marginalization over Qm- The 
solid, dashed, and dotted contours are for SNAP, LSST, and 
CFHLS surveys, respectively. Survey parameters are listed in 
Tableland a flat Universe is assumed. Fiducial models for 
both ACDM and a supergravity-inspired model are plotted 
as stars, and the shaded regions are the expected constraints 
from the SNAP Type la supernova measurement plus Qm 
prior (E. Linder, private communication). Unlike the weak 
lensing contours, the SN contours include the estimated ef- 
fects of the dominant systematic errors. 



cover « 200 deg^ of sky; a possible deep multicolor lens- 
ing survey covering ~ 10% of the sky with the pro- 
posed Large Synoptic Survey Telescope (LSST)^, and a 
400 deg^ weak lensing survey with the proposed Super- 
nova Acceleration Probe (SNAP)'^. The orbiting SNAP 
focal plane should obtain significantly higher source 
galaxy densities than the ground-based observatories, 
with consequently deeper median redshift. 

In evaluating the Fisher matrix, we assume that <7e — 
0.3 and that (3^ — 0.8. Cosmological parameter un- 
certainties will scales as (Je//3- We use Azi = 0.2 for 
the thickness of the lens slices. This choice has very 
little impact upon the parameter constraints because 
the first term dominates Equation (|22|l . e.g. changing 
to Az£ = 0.1 for the SNAP case alters the parameter 
constraints by < 1%. 

For each candidate survey we evaluate the 
cosmological-parameter covariance matrix Cp using 
Equation (|22|l . The resultant matrix is highly degen- 
erate in one direction, so we adopt a Gaussian prior 
distribution on Qm with a — 0.03. The constraints 
in the wq — Wa plane after applying this prior and 
marginalizing over are plotted in Figure |21 for the 
three surveys. Table ^ also lists the l-cr uncertainties 
in Wo and Wa (assuming in each case marginalization 
over all other parameters) and the correlation coefficient 

^ www. dnitelescope . org/dark_home . html 
^ snap.lbl.gov 



between these two parameters. 

For both the SNAP and LSST surveys, the constraints 
on the dark energy equation of state and its evolution 
are quite strong, comparable to or tighter than those of 
any proposed experiment of which we are aware. The 
l-cr uncertainties on Wa, even after marginalizing over 
all other parameters, are a^^ ~ 0.3 in both cases. In 
terms of the commonly used parameter w' = [dw/d ln(l + 
z)]z=i = Wa/'2 we have a^' ~ 0.15. 

3.5. Dependence of Constraints on Survey Parameters 

The Fisher uncertainties from the weak lensing survey 
will scale as /gky''^- Because the prior on Q.^ is indepen- 
dent of fskyi the marginalized error ellipses do not quite 
scale as /gky^^, but this is a valid approximation in the 
neighborhood of the canonical survey parameters. 

Increased survey depth leads to both greater source 
density n and greater median redshift Zmcd- Parame- 
ter uncertainties scale as (again not quite because 
of the prior), and hence unlike power-spectrum tomog- 
raphy, there is no breakpoint at which cosmic variance 
limits a survey of a given size. The additional depth is 
also a benefit: changing Zmod from 1.0 to 1.5 decreases 
the error on Wa by 23%, if the galaxy density and survey 
area are held fixed. In other words, the increased depth 
makes each galaxy about 1.5 times more valuable. 

The lower panel of Figuren]illustrates the contribution 
of different lens planes to the constraint on Wa for the 
SNAP survey. Most of the information comes from lens 
planes near z = 1, but if a deeper source distribution 
were available, lens planes at higher z would continue 
to add significant information. It would be particularly 
interesting to see how this behavior is altered with in- 
clusion of the cosmic microwave background as a source 
screen at z ~ 1000. 

Combining these various scalings, the overall con- 
straint on evolution of the dark energy equation of state 
will scale roughly as 

3.6. Comparison with Jain & Taylor 

Our parameter constraints can be compared with those 
obtained by Jain & Tavlor (2003) who used a simplified 
implementation of the cross-correlation approach. The 
measurement suggested by them was the tangential shear 
around foreground halos, identified using galaxy groups 
and clusters. They assumed that essentially all halos out 
to z = 1 with mass > 4 x IQ^^Mq could be identified 
this way. For a given lens slice they used source galaxy 
shapes on only ~ 10% of the sky. This approach cor- 
responds to taking a particularly simple construction of 
the template shear map. Other simplifications made by 
Jain & Taylor were that intrinsic ellipticity shot noise 
was taken to be the only source of error and only two 
bins in source redshift were used for each lens slice. We 
find that the method presented here can improve param- 
eter constraints, primarily due to the use of more than 
two redshift slices. However our results for tr^ and a^^ 
are close to those of Jain & Taylor because the size of 
the single-parameter marginalized errors is, for the large 
surveys considered, controlled by the prior for 
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Table 1. Future Weak Lensing Surveys and Dark Energy Constraints 



Survey 


Median 2 


Galaxy Density 


/sky 


wo Error 


Wa Error 


wo — Wa Correlation 






(arcmin~^) 








CFHLS 


1.0 


30 


0.005 


0.26 


1.04 


-0.92 


LSST 


1.0 


30 


0.1 


0.06 


0.30 


-0.80 


SNAP 


1.5 


100 


0.01 


0.10 


0.36 


-0.83 



4. SYSTEMATIC ERRORS 

The cross-correlation cosmography technique is essen- 
tially insensitive to several of the largest systematic error 
sources in power-spectrum tomography, namely residual 
PSF distortion and errors in calculation of the non-linear 
power spectrum. The technique does, however, place 
stringent demands on the distortion calibration and pho- 
tometric redshift accuracy of the weak lensing survey. 
This is illustrated in Figure ^ The top panel plots the 
fractional change in gis when we move from a ACDM 
universe to one with Wa — 0.2. According to the Fisher 
calculation, this should be detectable at la significance 
in the SNAP survey. 

4.1. Demands on Distortion Calibration 

Each line in the plot shows g^s vs Zg for one choice 
of zi . The "signal" in the cross-correlation cosmography 
method is the departure of these curves from horizontal 
lines. A pure vertical shift of any line will be degener- 
ate with a change in the fidelity /?£ of the template in 
lens shell £, hence the cosmological information is in the 
slope/curvature of these lines. We see that the mini- 
mum detectable cosmological signature corresponds to a 
change of « 2 parts in 10'^ of g, and hence an equiva- 
lent change in the measured distortion vs Zg. We may 
immediately conclude that exploitation of this technique 
to constrain Wa requires that the distortion calibration be 
constant to 1 part in 10^ over all measur ed redshifts. This 
will b e a significant technical challenge. iHirata fc Seliald 
()2003|) demonstrate that present algorithms have calibra- 
tion errors of 1-10%, so at least an or der-o f -magn itude 
improvement is required. [Hirata fc Seliald l|2003j) and 
iBernstein k, Jarvi^ l)2002() propose methods to improve 
calibration accuracy, but these have not yet been demon- 
strated on real data. Furthermore, most methodologies 
assume that the distortion is weak, i.e. they ignore in- 
duced changes to galaxy shape that are of order or 
higher. The non- linear effects of even "weak" distortions 
will have to be accounted for to second or even third or- 
der to reach the calibration accuracy of 10"'^. Regions 
where the lensing is strong will clearly have to be avoided. 

To estimate the demands of this requirement on the 
data reduction methodology, we note that the conver- 
sion of observed ellipticities to pre-seeing elliptiticities 
typically involves the factor (1 -t- rj/r^), where r is the 
angular radius of the (pre-seeing) galaxy and r^, is the 
PSF radius. A good weak lensing methodology makes 
use of galaxies as small as r « r^,; in this case, both r 
and r^. must be known to better than 1 part in 10^ in or- 
der to obtain a distortion calibration accurate to 1 part 
in 10~^. If r were independent of z, then we would have 
redshift-independent calibration errors, which do not af- 
fect the cross-correlation method. But fainter, more dis- 
tant galaxies are typically smaller in angular diameter, so 



errors in r^, will couple to the cosmography measurement. 

4.2. Demands on Redshift Determinations 

Our formalism can be adapted to deal with the ran- 
dom errors in photometric redshift measurements, but it 
will still be crucial to minimize biases in the photo-z esti- 
mates. The quantity gis scales slightly less than linearly 
with the redshifts Zg and zi. It is therefore clear that 
mis-estimates of the mean photometric redshift of a few 
parts in 10^ would overwhelm the signal of a = 0.2 
cosmology. Hence photometric redshifts must be accurate 
to Ri 10~^ in log(l -I- z). This again is at least one order 
of magnitude beyond the present state of the art. 

It is important to realize first that this is not the re- 
quirement on the accuracy of each measured photometric 
redshift, but rather a requirement on the bias of a collec- 
tion of photometric redshift estimates. Second, we are 
not required to use all the galaxies in the image. One 
would likely choose to exclude from the analysis galaxies 
whose colors make assignment of a photometric redshift 
particularly uncertain or ambiguous. 

4.3. Practical Issues 

If the required accuracy of distortion and redshift cal- 
ibration cannot be achieved, one could introduce addi- 
tional free parameters in the model to represent calibra- 
tion errors. These terms would be purely functions of 
Zs, and hence in theory distinguishable from the cosmo- 
logical signals, which couple zi and Zg through ggg. The 
constraints on dark energy would be degraded, to an ex- 
tent that is calculable with further Fisher analysis. 

We note that if the random error in a typical photo- 
metric redshift is az ~ 0.03(1 -I- z), then it will take a 
spectroscopic survey of ~ 10 galaxies in order to check 
that the mean error (bias) is below 10^*^(1 + z). Repeat- 
ing this bias check in bins of redshift from 0-3 would thus 
require a sample of 10^ or 10^ spectroscopic redshifts to 
compare with the photometric redshifts. Redshift sur- 
veys of this size to limiting magnitudes of 24 or 25 will 
be feasible: indeed the DEEP2 survey^ is already well on 
its way to its goal of 65,000 spectra to Iab < 24.5. 

Given the extreme demands that cross-correlation cos- 
mography will place on the calibration of both the lensing 
distortion and photo-z calibrations, there will be a sig- 
nificant advantage to space-based observations. Ground- 
based analyses must deal with a constantly-varying PSF 
and atmospheric transmission function which will make 
it more difficult to achieve these accurate calibrations. 
A space-based platform will have the additional advan- 
tages of thermal stability, a much sharper PSF, and the 
possibility of using near-IR filter bands to improve the 
photo-z accuracy. 

^ http://deep.berkeley.edu 



8 



5. CONCLUSION 

We have presented a form alism for implementing the 
idea of iJain fc Tavloi) l)2003D to constrain cosmology by 
tracing the dependence of induced distortion on back- 
ground redshift for a fixed foreground mass template. 
Cross-correlating the background distortion with a se- 
ries of foreground mass templates has significant advan- 
tages over power-spectrum tomography, particularly its 
immunity to spurious distortion signals, and the ability 
to use non-linear lensing structures without having to 
model accurately the non-linear evolution of matter in 
the Universe. We believe this formalism makes use of 
all the information available in all lens-source pairs in a 
nearly optimal fashion; we find that a survey with median 
redshift ~ 1.5 and 100 galaxies per square arcminute can 

— 1/2 

constrain the dark energy parameters to cr^(, w 0.01/gj,y 

and Gu]^ ~ 0.035/j,jjy^^ after application of a practical 
prior on Q.^- 

Realization of the full potential of the cross-correlation 
cosmography method will require that techniques for the 
calibration of lensing distortion and photometric redshift 
be improved by at least an order of magnitude from 
present state of the art. There are no known fundamental 
barriers to this, but it will not be easy. We discuss ways 
of fitting for parameters in the calibration from the data. 
To reduce possible biases in the photometric redshifts to 
an acceptable level would require spectroscopic redshifts 
for 10"* — 10^ galaxies over the redshift range used in the 
analysis. An orbiting observatory may be preferred for 
obtaining the precision required for our method due to 
its greater photometric and optical stability, and access 
to the near infrared. 

The potential accuracy of cross-correlation cosmogra- 
phy on the equation-of-state time variation compares well 
to the expected precision of power-spectrum tomography. 
A formalism for full utilisation of the power-spectrum to- 
mography information has not yet been published, and 
none of the few published investigations of the effects 
of time-dependent dark energy use Equation (|23|l for the 
equation of state. The assumed priors and confidence lev- 
els of published estimates also vary gr eatly, so only crude 
comparisons can be made. |h3 l!2002br estimates an error 

of ~ O.OS/gi^y on w' from tomography confined to the 
linear regime. The plots of iBenabed fc van Waerbek^ 
l|2003(l suggest w 0.03/^^^^^ from a non-tomographic 
analysis deep into the non-linear regime. It should ul- 
timately be possible to use all the information encoded 
in the weak lensing — cosmography, growth function, and 
cross-correlation with foreground structures — to provide 
constraints stronger than those we forecast here. By can- 
celling the growth information, however, we eliminate the 
systematic errors that might arise from mis-calculation 
of the theoretical non-linear spectrum, and from addi- 
tive contamination of the distortion field by systematic 
effects. The penalty for cancelling the growth factor is 
that the signal is partially cancelled as well, leaving the 
cross-correlation method more susceptible to multiplica- 
tive, redshift- dependent errors in the distortion field that 
might arise from PSF effects. 

Further improvements to the cross-correlation cosmog- 
raphy methodology are worth investigation. The cos- 
mic microwave background can serve as an additional 



source plane at Zs = 1000, and its shear pattern can 
be cross-correlated with all the lens planes to provide 
a S/N improvement and greater redshift leverage. One 
could also make use of magnification information as well 
as shear in the background galaxy samples, which would 
improve the S/N and serve as a useful cross-check {e.g. 
IQain (2002)). Finally, in this study we have not devel- 
oped a detailed method to obtain template distortion 
maps from the foreground galaxy distribution. A de- 
tailed study that includes the redshift dependence of the 
fidelity of the template would be useful. 
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APPENDIX 

COVARIANCE MATRIX OF R^fs 

At several points in the analysis we require the covariances of the slice-to-slice distortion correlations Ra/j = (d^ • d^) . 
Under the assumption that d^ is completely uncorrelated with both d^ and d^, it is clear that the only non- vanishing 
elements of the covariance matrix are of the form Va.r{Raa), Var(i?Q,^), and Gav{RapRi3a)- We first calculate Var(i?ct^) 
for a ^ (3: 

{Rai3Rap) = -^ J J d^rdV (^[d^{r)-dp{r)][d^{r')-dp{r')]) (Al) 
- ^ / Ka+ {r)((}+ {r) + ^c- {r)^p- (r)] (A2) 
^ ^ '^'^P„(^)P^W. (A3) 



n J (27r) 

fl is the solid angle of the survey, and -Pq(^) and P0{£) are the power spectra of the lensing convergence corresponding 
to the distortion fields do, and d_a, respectively, at spherical harmonic £. We have made use of the shear correlation 
functions ^± as defined for example in Schneide r, van Waerbeke. fc Mellie^ l)2002fl . who also give these correlation 
functions in terms of the convergence power spectrum. We may also put the distortion variance in terms of the 
convergence power spectrum: 

(4) = ^'^^''(^)- 

We can now express the covariance straightforwardly as 

{R.0Rap)^{dl){4)^ (A5) 

"^^^ 2 jdHP^{e)jdHPp{e)- ^ ' 

The quantity f2ct^ gives the solid angle over which the two distortion fields maintain mutual coherence. We will make 

the approximation that the template field d^ has a power spectrum with the same shape as the actual field d^ , so that 
flap has the same value regardless of whether we are correlating the templates or the real distortion, e.g.Q.pa — ^af3- 
If the lens shell thickness Ax satisfies Ax «C Xf; then we have a simple relation between the convergence power 
spectrum and the spectrum of mass fluctuations in three dimensions: 

m) - !Tfh p^<^ i -^Xp] Ax (A7) 



4a2(x/3)c'* VX/3 

Cav{RapRpa) is also non-vanishing; it differs from Var{i?c(^) in that cross-correlations between d^ and d^ are 
required. If the template is well correlated with the actual mass then we expect Equation IjASp to describe the 
covariance up to a factor near unity, i.e. R^fs will be highly correlated with Rf3a. We do not need this expression in 
detail. 

The final non-vanishing term is (again ignoring the distinction between template and actual distortion) 

Var(i?„o) = l J d\{[dl{0)-(d')] [dl{r)-{d')]), (A8) 

which clearly involves the four-point correlation functions of the shear field, and hence depends upon the degree of 
non-Gaussianity in the mass distribution. Equation (jA5|l will be correct up to a factor of order unity, which we could 
absorb into the definition of ^laa- This term does not significantly affect the Fisher matrix. 



